arXiv:1509.05804vl [cond-mat.mes-hall] 18 Sep 2015 


Theory and experiment on cavity magnon polariton in the ID confignration 

B. M. Yao/’2-E] Y. S. Gui/ Y. Xiao,^-^ H. Guo,'^ X. S. Chen,^ W. Lu,^ G. L. Chieii,^ and G.-M. Hu^ 

^Department of Physics and Astronomy, University of Manitoba, Winnipeg, Canada R3T 2N2 
^National Laboratory for Infrared Physics, Chinese Academy of Sciences, Shanghai 200083, People's Republic of China 
^College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China 
‘^Department of Physics, McCill University, Montreal, Quebec, H3A 2T8, Canada 
^Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA 

(Dated: September 22, 2015) 

We have theoretically and experimentally investigated the dispersion of the cavity-magnon- 
polariton (CMP) in a ID configuration, created by inserting a low damping magnetic insulator into 
a high-quality ID microwave cavity. By simplifying the full-wave simulation based on the transfer 
matrix approach in the long wavelength limit, an analytic approximation of the CMP dispersion 
has been obtained. The resultant coupling strength of the CMP shows different dependence on the 
sample thickness as well as the permittivity of the sample, determined by the parity of the cavity 
modes. These scaling effects of the cavity and material parameters are confirmed by experimental 
data. Our work provide a detailed understanding of the ID CMP, which could help to engineer 
coupled magnon-photon system. 

PACS numbers: 85.75.-d, 71.36.+C, 76.50.-|-g 


I. INTRODUCTION 

Coupling between electrodynamics and magnetization 
dynamics is a subject of long-standing intereslii^—. At 
the resonant condition where the microwave frequency 
approaches the eigenfrequency of a magnon, two pieces 
of key physics stand out from the resonant magnon- 
photon coupling: (i) mode hybridization takes place 
which leads to the formation of a magnon polaritort^, 
(ii) damping correlation happens which causes radiation 
dampingiii^i^d^i^^ . 

Very recently, interest in the physics of magnon-photon 
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TABLE I. Summary of different theoretical approaches de¬ 
veloped receirtly for studying the physics of magnon-photon 
coupling. 


coupling has grown significantly due to the advancement 
of microwave cavity and spintronic techniques. Experi¬ 
mental progress has shown that the strong coupling be¬ 
tween magnon and photon can be easily achieved by 
inserting a low damping magnetic material into a high 
quality microwave cavity^^— . Such a cavity spintronic 
methodic creates new avenues for studying and utiliz¬ 
ing both spintronics and quantum information, as high¬ 
lighted by the recent development of a cavity spin pump¬ 
ing method for coherently manipulating spin currents us¬ 
ing magnon-photon couplingi^ and the invention of quan¬ 
tum transducers that coherently link diverse quantum 
systems^SiSi, respectively. 


Approaches 

Classical theory 

Quantum theory 

Intuitive 

coupled harmonic oscillators [Bai et al., Ref. 

coupled harmonic oscillators [Zhang et al.. Ref. ^ 

Concise 

electrodynamic phase correlation [Bai et al., Ref. 

quantum entangled states [Soykal et al.. Ref. ^3] 

Detailed 

ID scattering theory [Cao et al., Ref. 2^ 

transfer matrix method [this work] 



Theoretically, as summarized in the Table I, different approaches have been developed for studying the reso- 













2 


nant magnon-photon coupling. The simplest approach 
roots on the physics of coupled harmonic oscillators, 
which provides an intuitive picture of the key physics 
of mode anticrossing and damping evolution. This ap¬ 
proach can be formulated by using either a classical or 
a quantum mechanical model, as detailed in the Supple¬ 
mentary Materials of Ref. and[l^, respectively. Within 
such an approximation of harmonic oscillators, there is 
no distinction between the physics of quantum and clas¬ 
sical coupling. 

Such a distinction appears in the refined general ap¬ 
proaches which specifically identify different origins of the 
magnon-photon coupling in the quantum and classical 
regimes. In the fully quantum-mechanical treatment^^, 
it was found that the quantum-coherent magnet-photon 
coupling arose from the entangled quantum states of 
spin orientation and photon number. On the other 
hand, in a concise but general classical treatment in¬ 
dependent of the sample and waveguide geometry^, it 
was found that classical-coherent magnet-photon cou¬ 
pling appears due to the mutual electrodynamic cou¬ 
pling of the macroscopic microwave field and the macro¬ 
scopic magnetization, which follows the classical phase 
correlations determined by Faraday’s law and Ampere’s 
law (i.e. by Maxwell’s equations). So far, it remains 
a theoretical challenge of clearly define the border be¬ 
tween the quantum and the classical coupling regimes, 
which is needed for understanding the transition from 
quantum- to classical-coherent magnet-photon coupling. 
However, in light of the general classical modeli^ which 
quantitatively captures both key features of mode hy¬ 
bridization and radiation damping measured in recent 
experiments^^— ^ it appears that many of these exper¬ 
iments probe the property of cavity magnon polari- 
tons (CMPs) which are formed by the classical-coherent 
magnet-photon coupling. 

As known from early theoretical study of the propaga¬ 
tion of electromagnetic waves in a microwave waveguide 
partially filled by a magnetic slab^, the details of the 
CMP tediously depend on the parameters of practical 
interest, such as the shape and size of both the waveg¬ 
uide and the sample. In general, it requires numerical 
solution of the coupled Maxwell’s equations and Landau- 
Lifshitz-Gilbert equation. Only in special cases, simple 
solutions can be obtained. One of such special examples 
is the ID cavity inserted with a magnetic slab, which 
was recently theoretically analyzed by using the ID scat¬ 
tering theory^!. To verify the specific theoretical results 
obtained from such a detailed approach, an experimental 
method of studying the classical-coherent magnet-photon 
coupling in the ID configuration needs to be developed. 

In this paper, we develop such a ID method experimen¬ 
tally. Furthermore, we simplify the ID scattering the¬ 
ory by using the straightforward transfer matrix method, 
so that our experimental and theoretical method can be 
easily used in a combined way for studying the detailed 
physics of the CMP in the ID configuration. One of our 
focus here is to explain how the coupling between pho¬ 


tons and magnons depends on the thickness of the sample 
as well as the parity of cavity mode for a ID cavity. 

This paper is split into two main sections: theoretical 
model and experimental results. In the theoretical model 
section, we first provide a brief description of the ID cav¬ 
ity used in this work and use the transfer matrix approach 
to derive the cavity resonance. Then the dispersion of the 
CMP due to the hybridization of the cavity resonance 
and the ferromagnetic resonance (FMR) in a bulk yt¬ 
trium iron garnet (YIC, Y 3 Fe 50 i 2 ) is derived by solving 
the coupled Maxwell’s equations and Landau- Lifshitz- 
Cilbert equation using the transfer matrix method. To 
highlight the coupling strength in a CMP, the full-wave 
simulation is simplified in the limit of long wavelength, 
which clearly shows its strong dependence on the cavity 
length, sample thickness, sample permittivity as well as 
the parity of the cavity modes. Finally we present ex¬ 
perimental results to confirm our numerical results and 
analytical model. 


II. THEORY OF CAVITY MAGNON 

POLARITON IN THE ID CONFIGURATION 

For a quantitative understanding of the CMP disper¬ 
sion, we developed a model to solve the case for a ID 
microwave cavity. In general it is difficult to provide the 
exact analytical solution of modes excited in a ferrite 
loaded cavity under a magnetic bias, where the propa¬ 
gating wave encounters spin dynamics in the loaded fer¬ 
rite sample^ i^^:^^ . However, in the special case of the one 
dimensional cavity that is used in this work the prob¬ 
lem of wave propagation can be analytically solved— 
To simplify the complicated mathematical process, the 
scattering matrix treatment and transfer matrix, being 
derived from general principles such as conversion of en¬ 
ergy, are used to deduce the exact solution of the CMP. 
This full-wave solution provides all detailed information 
of the coupled magnon-photon system without any ad¬ 
justable parameters. In addition, a simplified analytical 
solution is also deduced in the long wavelength limit and 
validated by the full-wave calculation, allowing the pre¬ 
diction of the coupling strength between cavity modes 
and magnon modes in a straightforward way. 


A. Microwave cavity based on waveguide assembly 

As shown in Fig. [Ija) the microwave cavity used 
in this work is a Fabry-Perot-like resonator based on 
the Ku-band (12 CHz- 18 CHz) assembled waveg¬ 
uide apparatu a^^’^^'^^ , where circular waveguides are 
connected through circular-rectangular transitions to 
coaxial-rectangular adapters and two transitions are ro¬ 
tated by an angle 6 . The microwave propagation in such 
an apparatus can be characterized by the S-parameter 
and measured by a vector network anlayzer. 

The fundamental mode of the rectangular waveguide 
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(a) 



(b) Port A Port L Port R port B 



FIG. 1. (Color online) (a) Schematic of the waveguide assem¬ 
bly cavity. Two circular-rectangular transitions are twisted 
by an angle 9 with respect to each other and are connected 
to a vector network analyzer through SMA to rectangular 
waveguide adapters, (b) Wave propagation between port A 
and B in (a) described by a cascade connection of two-port 
networks, (c) Schematic diagram of the wave propagation for 
incident wave from port 1 and transmitting wave toward port 
2. For the coupling experiments, the static magnetic bias H^xt 
is applied in the x — y plane, perpendicular to the direction 
of wave propagation (along z direction). 


is TEIO and the polarization of the microwave electric 
held is along x and x' (the short axis of the rectan¬ 
gle in Fig. UKa)) for port A and B, respectively. For 
the circular waveguide the fundamental mode is TFll, 
supporting two degenerate orthogonal waves, named x— 
and y—polarized waves for port A {x' — and y'—polarized 
waves for port B). Therefore, a circular-rectangular tran¬ 
sition is designed to smoothly transform the TEIO mode 
of its rectangular port (port 1) to the co-polarized TEll 
mode of its circular port (port A) without any reflection. 
However, propagating from port A towards port 1 only 
x-polarized waves be can transmitted through while y— 
polarized waves are rehected with a rehection coefficient 
of R (very close to 1) and a phase shift of (j)y, which can 
be determined experimentally. Similar effects appear for 
x'- and y'— polarized waves propagating form port B 
toward port 2. Because of this mirror-like effect for y— 
and y' — polarized waves, a Fabry-Perot-liked resonator is 
built and a standing wave is generated in the long dimen¬ 
sion (along the z direction) of the circular waveguide at a 
particular frequency, which can be deduced by Maxwell’s 
equations as detailed in the discussion below. 


B. Analytical model of the resonant cavity 


tive analysis of the coupling between magnon and cavity 
modes, we start with the transfer matrix (the ABCD 
matrix in microwave terminology) for a dielectric loaded 
cavity. As shown in Fig. [TJb), Ig is the thickness of 
the sample loaded waveguide and I is the length of air- 
filled waveguide. Hence the whole length of the circular 
waveguide is W=2l+ls. Following the convention, 

the transfer matrix for a two-port network, for example, 
the sample-loaded waveguide (yellow area between port 
L and R in FiglTJb)), can be defined as^S 

feA _ ( cos(44) -iZsS\n{kglg)\ f e/j\ _ f cr 

\-iZA sin{hQ cos{kglg) J ArJ 

( 1 ) 

where the subscript L{R) indicates the left(right) inter¬ 
face of the sample. Equation © links the electric held 
(e) and magnetic held (h) of microwave propagation in 
a dielectric medium using the ABCD matrix M^, where 
the matrix elements relate only on the material param¬ 
eters of wave number (fc, = (lu/ c)A^sl^s — (wc/w)^), the 
characteristic impedance {Zg = ujyofjg/kg) and the length 
of transmission line (4). Here w, tq., £s(ys), and £o(mo) 
are the frequency of the operating microwave, the cut-off 
frequency of the circular waveguide, the relative com¬ 
plex permittivity (permeability), and vacuum permittiv¬ 
ity (permeability), respectively, c = 1/^eoMo is the vac¬ 
uum speed of light. For the air-hlled region in waveguide, 
the wave number and characteristic impedance are fco and 
Zq using £g = yg = \. 

We note that the hnite lateral size of the waveguide re¬ 
sults in the cut-off wavelength Ac = 27rc/wc; for the TEll 
mode in a circular waveguide Ac is 3.41 times the radius 
(8.05 mm for standard Ku-band circular waveguides)^, 
corresponding to a;c/27r = 10.85 GHz. As a consequence, 
the velocity of microwave propagation along the waveg¬ 
uide is less than its velocity through free space (speed 
of light), which results in a reduction of the ratio of the 
effective sample thickness with respect to the effective 
length of the cavity and hence significantly affects the 
dispersion of the CMP. 

The convenience of the transfer matrix relies on the 
fact that the ABCD matrix of the cascade connection of 
the two networks is equal to the product of the ABCD 
matrices representing the individual two-ports. As such 
in Fig. [ijb), the transfer matrix between ports A and 
B is Mab = Ml Mg M 2 and Mi and M 2 are the air-filled 
circular waveguide located at the left and right side of 
the sample, respectively. However, the ABCD transfer 
matrix relates the total electric and magnetic field at the 
ports; to characterize the fields incident on the ports to 
those reflected from the ports the scattering matrix (also 
referred to as S'-parameter) should be involved, defined 
as 


(e-A _ (Saa Sab\ (e+\ 
[e-J - [Sba SbbJ [e+J 


( 2 ) 


To provide the analytical solution for the wave propa¬ 
gation in the microwave cavity and hence the quantita- 


where the superscript indicate the waves enter¬ 

ing/exiting the port. The elements of the scattering ma- 
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trix can be easily obtained from the ABCD transfer ma¬ 
trix according to Eq. (IA7I) in Appendix The great 
advantage of using the scattering matrix is that the S- 
parameter can be measured directly with a vector net¬ 
work anlayzer. 

As the experimentally accessible parameter is the 


transmission coefficient S 21 between ports 1 and 2 of the 
assembled waveguide cavity, a wave entering or exiting 
the ports A or B must be considered by including the 
boundary condition. A detailed mathematical process 
can be found in Appendix!^ leading to the transmission 
coefficient S 21 from port 1 to port 2 


cos(0)5'ab[—1 + SabSba + b^’^'^RSbb — RSaa{—^ + RSbb)] 
— 1 -I- cos^{9 )R'^SabSbA + e^'^yRSsB — e’-'f’yRS aa{,— ^ + e^^^yRSBB) 


This equation is the key result of the full-wave solu¬ 
tion for calculating the property of CMP in the ID wave 
propagation configuration. Here R and (jiy are the re¬ 
flection coefficient and phase shift for the non-supported 
polarized microwave for the rectangular waveguide, re¬ 
spectively. For an air-filled waveguide (empty waveguide, 
Es = = 1), it is found that Saa = Sbb = 0 and 

Sab = Sba = and hence 

1 p2 2i(2fcoJ-|-fco4-l-0H'| 

U21 ^ (,Qg2(0)g2i(2fcoi-4fco4+0D ' 

. 

The waveguide assembly acts as a transmission 
through waveguide with 5'2i = at cos(0) = 1, 

while it would block transmission with S '21 = 0 at 
cos{9) = 0. Except for these extreme cases, the spec¬ 
tra of S '21 appears as a set of cavity modes occurring 
at a cavity resonant frequency of locr with a mini¬ 
mum transmission of |S' 2 i(wcb)| = cos(0)(l — RS)/[l — 
i?^cos^(0)] ~ 0, where kcRi‘21 + 4) + with 

kcR = i^CR/c)y/l - {lUc/uicr)'^, q being an integer. 

In our experiments, the angle of 9 is rotated to be 45° 
unless otherwise specified. As expected a set of reso¬ 
nances is observed for the air-filled waveguide as well as 
the Teflon-loaded waveguide shown in Fig. [2](a). These 
cavity resonances can be classified into two groups: e- 
modes (corresponding to even number q) being anti¬ 
nodes of the microwave electric field and h-modes (cor¬ 
responding to odd number q) being nodes of the mi¬ 
crowave electric field at the sample position. In the ab¬ 
sence of any external applied static magnetic bias, i.e., 
Hext = 0, the loaded dielectric sample only causes the 
red-shift of the e-modes but has negligible effect for the 
h-modes; this effect can be clearly seen in Fig. [UJa) when 
a 1.6 mm thick Teflon disk is loaded in the cavity. Us¬ 
ing the full-wave simulation base on Eq. the pre¬ 
cise values of R and (j)y were determined for the air-filled 
waveguide: for example, i?=0.995 and 4>y = —36.1° for 
the e-mode occurred at 13.199 GHz and i?=0.995 and 
cjiy = —6.5° for the h-mode occurred at 12.159 GHz. 

In order to distinguish the different dielectric depen¬ 
dence of e- and h-modes on the loaded material, we have 
simplified Eq. m near the cavity modes by Taylor ex¬ 
pansion. In the limit of long wavelength, to 11(^41 ^ 1, 
and assuming R = I for simplicity, one can obtain the 


first order approximation given by 




1 - 


1 


S,.=-^e-^*y 


1 - 


1 - i[6l{ko - kcR) + 34(^ - ho)] 
(e-mode) 

1 


1 - i[6/(fco - kcR) + 34(/4 - l)fco]. 

(h-mode) 

(5) 


The prefactor of S '21 is determined to be 2-\/2/3 since 
far from the cavity modes IS 21 I tends to be 21 cos(h)|/[l-|- 
cos^(h)] according to Eq. (g]). To clearly see the resonant 
feature of the cavity modes, we rewrote the above expres- 



FIG. 2. (Color online) (a) Measured amplitude of S 21 spectra 
for the air-filled cavity (solid line) and Teflon-loaded cavity 
(dotted line) at 0 = 45° and Hext = 0 with 1U=85 mm and 
Zs = 1.6 mm.(b) Calculated IS 21 I spectra according to Eq. (j6]) 
for air-filled cavity (Black line) and Teflon-loaded cavity (red 
line) as compared to the measured data (symbols). 
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sions as 


C. Full-wave solution of the CMP dispersion 


^21 

iS '21 ■ 


2^2 


M,-. 




i{uJ — UlCR + Su/cr) — ^^CR 

(e-mode) 

I_ Aa;cj; _ 

i{u) — LUCR + SlUcr) — ^^CR, 

(/i-mode) 


( 6 ) 


with 


Si^CR = 


^Jsj \UJCR 




\^crJ 


, (e-mode) 
(/i-mode) 

(7) 


Here the linewidth of the cavity mode is Alocr = 
(c/60yrr {uJc/ojcrY- Physically, the linewidth of the 
resonator is related to the energy loss of the cavity loaded 
in the circuit, which can be estimated by the lifetime 
of the photon trapped in the cavity. Due to the fi¬ 
nite cross section of waveguide, the velocity of the mi¬ 
crowave photon at resonances is C\/l — (wc/wcfi)^ along 
the waveguide, and hence the travel time is about = 
2llcy/l^ {oJc/ijJcrY ■ Depending on the polarization, 
half of the photons could be lost at the interface of either 
port A or B in our proposed ID cavity; thus, the average 

+ 00 

lifetime of photon is about n = Tt = 2 Tt, which 

n—1 

is the same order as I/Awch = 3rt estimated from the 
transfer matrix approach. 


Equations and © clearly indicate that the h- 
mode is not sensitive to the dielectric material (usu¬ 
ally ^ = 1 ) loaded at the center of the cavity where 
the node of the microwave electric field is present. Us¬ 
ing the values of £« = 2.0 (at 13 GHz provided by the 
company), wcfl/27r=13.199 GHz (12.159 GHz) for the 
e-mode (/i-mode), and 4 = 1-6 mm, the measured S- 
parameter for the dielectric loaded cavity can be well 
reproduced as shown in Fig. [UJb) based on the calcu¬ 
lation according to Eq. ([5)). The red-shift of the e- 
mode shows excellent agreement between measurement 
and calculation with Sujcr/“271=0.25 GHz. The calculated 
linewidth of AujcrI27t = {c/6l)^l — {ujc/locrY = 0.11 
GHz is slightly smaller than the measured linewidth of 
Aujcr/27t = 0.14 GHz, which may be due to the contri¬ 
bution from intrinsic damping caused by factors such as 
conductor loss of the cavity, not included in our model. 
The intrinsic quality factor of cavity mode which 

denotes the loss in the resonator itself can be estimated 
as ~ 10 ^ approximately according to the ratio between 
cavity resonance frequency and the bandwidth at -|-3dB 
with respect to the minimal transmission^^. 


Far from the FMR frequency lofmr of the ferrite sam¬ 
ple, the spin-precession frequency is either much larger 
or much smaller than the frequency (w) of the driv¬ 
ing microwave held. Therefore, the interaction between 
the electromagnetic held and spin dynamics is negligible. 
As a consequence, the YIG sample acts as a dielectric 
medium, only causing the red-shift of the e-mode while 
causing no hrst order effect for the h-mode as discussed 
in Sec. IHBI In the vicinity of the FMR frequency, the 
dynamic response of the spin precession m driven by the 
microwave magnetic held h is governed by the Landau- 
Lifshitz-Gilbert equation^. Meanwhile, the wave propa¬ 
gation in the ferrite should also satisfy Maxwell’s equa¬ 
tion. 

When UJPMR approaches lucr, the electrodynamics and 
magnetization dynamics are strongly coupled. Gon- 
ventionally, one can solve the coupled Landau-Lifshitz- 
Gilbert equation and Maxwell’s equation as detailed in 
the discussion in Appendix[B]by introducing the effective 
permeability fieff given by 


Me// — 


i^+XL^-X^T 

1 + XL 


( 8 ) 


where xl and xt are the diagonal and non-diagonal ele¬ 
ments of the Polder tensor as dehned in Eq. (B5) in Ap¬ 
pendix [B) The static magnetic bias Hext is applied along 
the y'-(;jii-ection, perpendicular to the direction of wave 
propagation (along z direction) in a coordinate shown in 
Fig. [H 

When luj/ujpMR — 1| <C 1, the resonant feature of ^eff 
can be expressed as 


.. ^ 1_ (cug -I- 

2 (w — WFMfi) + -I-a;M/2) ’ 

wit h cug = xHex t, OJM = 7 M 0 and ujfmr = 

1\/Hext{Hext + Mq). Here 7 = Hoe/mg is the effective 
gyromagnetic ratio of an electron (charge e and mass 
TOe), a is the Gilbert damping parameter, and Mq is the 
saturation magnetization. 

Replacing in Eq. (IA5I) with the effective permeabil¬ 
ity Heff, and substituting the resultant wave number kg 
and the characteristic impedance Zg for the YIG sample 
into the transfer matrix, the S-parameters {Saa, Sab, 
Sba and Sbb) between port A and B are deduced as 
detailed in the discussion in Sec. HI Then the transmis¬ 
sion coefficient near the FMR can be calculated by the 
full-wave solution according to Eq. ©■ With param¬ 
eters of Wc=0, (/>y=0, i?=l, fioMo = 175 mT, 7 = 176 
fiQ GHz/T, a = 1.25 x 10“^, and W = 46 mm. Figure 
[3Kb) plots the calculated 5 ' 2 i spectra as a function of mi¬ 
crowave frequency and magnetic field near the h-mode 
(g=3) occurred at uj/27t =9.77 GHz. Without going into 
the detailed behaviour of the CMP, we focus on the cou¬ 
pling strength {g) of the CMP, which is determined by 
the CMP gap occurring at ojfmr = <^CR as shown in 
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(a) 



ISj.UdB) 

(C) 


(b) 




H„H„, (mT) 



FIG. 3. (Color online) (a) S 21 spectra at lofmr = ujcr=9-77 
GHz, where the CMP gap was used to determine the cou¬ 
pling strength g. (b) Dispersion of the CMP as a function 
of microwave frequency and magnetic field at ls=5 g.m. The 
dashed line indicates the condition of ufmr = ucr- (c) Cou¬ 
pling strength with q=3 and 4 are calculated based on the 
scattering theory (open symbols, from Ref. [ 2 ^ and trans¬ 
fer matrix theory (solid lines, this work) with parameters of 
tJc=0, hqMo = 175 mT, 7 = 176 /ro GHz/T and IF = 46 mm. 


Fig. IHa), to highlight the general feature of the coupling 
between the magnon and the photon in a ID cavity. 

Under such an approximation of zero cut-off frequency 
(ojc = 0), we compare the results from different ap¬ 
proaches: the scattering theorj*^! and transfer matrix 
theory. The calculated cavity modes corresponding to 
q = 3 and g = 4 for the empty cavity occur at 9.77 GHz 
(h-mode) and 13.03 GHz(e-mode), respectively, which 
are very close to the results obtained from the scattering 
theory. The coupling strength (solid lines) for both h- 
mode and e-mode are plotted as a function of IJ shown 
in Fig. He). Despite significant differences in the mi¬ 
crowave feed for the cavity used in the scattering the¬ 
ory (open symbols) and transfer matrix theory (solid 
lines) to model the GMP, the final result of the coupling 
strength is consistent, indicating the universality of the 
CMP, which is independent of the specific design of the 
ID cavity. 

In a practical implementation, the impact of the cut¬ 
off frequency should be taken into account because of 
the hnite lateral size of the ID cavity. To demonstrate 
this effect, the coupling strength were calculated at 07=0 
and 07=10.85 GHz (determined by the size of the cir¬ 



FIG. 4. (Color online) The coupling strength for neighbour¬ 
ing /i-mode (circles) and e-mode (squares) in a ID cavity 
(with 1U=85 mm) for oic=0 (open symbols) and oic=10.85 
GHz (solid symbols). For 10^=0, /i-mode (e-mode) occurs at 
12.34 GHz(14.11 GHz) at = 0 and R = 0.995, while for 
o)c=10.85 GHz, /i-mode (e-mode) occurs at 12.28 GHz(13.20 
GHz) at (f)y = —36.5° and R = 0.995. The dashed lines are 
guided by oc and oc z/^ relations for the h- and e-modes, 
respectively, in the limit of the long wavelength approxima¬ 
tion. 


cular waveguide in our cavity) by varying Ig more than 
three orders of magnitude from 1 fj,m to 1.5 mm. For 
the Zi-mode, the coupling strength has about a factor of 
2 reduction in the entire calculated Ig range when using 
a;c=10.85 GHz as shown in Fig. 0] In contrast, the cou¬ 
pling strength for e-mode is not sensitive to 07 in the 
limit of long wavelength Ig < 0.2 mm (indicated by the 
grey area in Fig. H), in which {./eIuicr/ c)lg < 0.2 with 
eg Ril5 for the YIG sampl o^^i^^ . 

Figure |4] also shows that the coupling strength for the 
Zi-mode and e-mode follows different power dependence 
of Ig guided by the dotted lines in the limit of long wave¬ 
length. Beyond this range the red-shift of the e-mode 
causes its approach to the neighboring Zi-mode and the 
hybridization of e- mode and Zi-mode due to the dielec¬ 
tric coupling in the cavity may complicate the determi¬ 
nation of the intrinsic coupling strength of the magnon- 
photon system. In addition, the cavity damping param¬ 
eters A.UJCR of both Zi-mode and e-mode are strongly de¬ 
pendent on ujc (not shown). The quantification of cavity 
damping is very important to characterize the phenom¬ 
ena associated with distinct coupling range for the cou¬ 
pled magnon-photon system. 


D. Analytic solution of the coupling strength 

To verify the different scaling effects for h- and e- 
modes and hence provide an analytical solution of the 
coupling strength between magnons and photons in a ID 
cavity, we have simplified the dispersion of the CMP in 
this section. To do so. Equation (H is substituted into 
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Eq. (O and the Taylor expansion of 4 is used to de¬ 
termine the dominant power dependence of the coupling 
strength. 


1. Coupling strength g for h-modes 


In the limit of long wavelength, the first order approx¬ 
imation is sufficient to model the dispersion of the CMP 
associated with /i-modes. In this case, the transmission 
coefficient near the FMR is reduced to 


2\/2 

5^21 =- 




i{uj - ujcr) - Aujcr + ■ 


( 10 ) 


with 

Awcfl = {c/6l)^/^ — {uJc/uJcrY: 

^^FMR = O!{l0h + UJmI"^), 

9^ = (4/40[l — (wc/wcfl)^](wff-I- wm)wm, (11) 

where Aojcr is the damping coefficient for the cavity res¬ 
onance, Aujfmr is the damping coefficient for FMR and 
g is the coupling strength between them. In Eq. (HU, the 
coupling strength g for /i-mode is linearly proportional to 
and independent of the permittivity of the 
sample. It was also found that the cut-off frequency Wc 
causes significant reduction of the coupling strength g of 
the /i-mode as well the damping of the cavity. 

To validate the analytical expression Eq. (HU for the 
/i-mode, we compare it with the results using full-wave 
simulation based on Eq. ©• Parameters for our YIG 
sample are gyromagnetic ratio of 7 = 169 /tq GHz/T, the 
saturation magnetization of goMo =0.155 T and damp¬ 
ing factor a = 1.25 x 10“^ obtained from independent 
measurements, and Eg = 15 is from Ref. [13 and HE As 
shown in Fig. [Ua), the coupling of the cavity and FMR 
modes results in a CMP gap i=2g) at ujcr = ^fmr, 
where the magnon-like CMP and the photon-like CMP 
have an equal amplitude in S 21 spectra. 

For a systematic comparison, the full-wave simulation 
is performed as a function of 4, I and Eg at the con¬ 
ditions of LucR = wfmr=12.16 GHz (corresponding to 
Mol^ext=380 mT) and the resultant coupling strength g 
is plotted as solid symbols in Figs. lUb-d), respectively, 
where IF = 2/ -|- lg=85 mm and es=15 in Fig. EKb), 
4=0.03 mm and es=15 in Fig. EKc), and 4=0.03 mm and 
IF = 85 mm in Fig. [51(d). In order to satisfy the condi¬ 
tion of a;cij=12.16 GHz the odd integer q and (fy should 
be adjusted when varying I during the full-wave simu¬ 
lation. However, both should not affect the calculated 
coupling strength. The coupling strength g calculated 
from the analytical expression of Eq. HU is also plotted 
as a solid line in Fig. |SKb-d). The agreement between our 
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FIG. 5. (Color online) (a) Simulated S 21 spectra of /i—mode 
for the YIG loaded cavity at goHiixt=380 mT (red solid lines) 
as compared to the cavity mode for the empty cavity (black 
dotted lines). The coupling strength <7 as a function of (b) 
Zy^ at IF = 85 mm and Es = 15, (c) Z“b2 — q Q 3 
and Es = 15 and (d) Ss at IF = 85 mm and h = 0.03 mm. 
Solid lines are the analytical results based on Eq. HU and 
solid circles are extracted from the full-wave simulation based 
on Eq. (|3|. 


full-wave simulation and analytical approximation is re¬ 
markable. Therefore, Equation HU provides a rigorous 
relationship between the sample parameters, the geomet¬ 
ric parameters of the resonant cavity and the coupling 
strength, which should benefit not only the study of the 
GMP but also the design of related microwave devices. 


2. Coupling strength g' for e-modes 

For the first order effect, the loaded YIG sample will 
only result in the red-shift of e-modes, which is deter¬ 
mined by the electromagnetic property of the sample, 

expressed as Sujcr = ^w[(ei — 1 ) + (1 -)(——)^] in 

2Z pig LOCR 

Eq. d?]). This effect is illustrated in Fig. UKa), where the 
full-wave simulation indicates that the cavity resonance 
for the e—mode originally at 13.20 GHz shifts to 13.00 
GHz when inserting the YIG sample in the absence of the 
external magnetic held. More interestingly, the full-wave 
simulation clearly reveals the resultant CMP gap due to 
the hybridization of the cavity modes and FMR modes at 
= 412 mT (corresponding to ujfmr = ojcr - Sujcr)- 
To analytically explain this effect, the effect due to higher 
order Ig terms should be taken into account. Differ¬ 
ent from the analytical approximation for Zi-modes, one 
should expand Eq. m according to both w — ujcr and Ig , 
because in this case uj — ujcr ~ Sujcr, i.e., oc 4- Expand¬ 
ing up to third order in both uj — uicR and 4 and replacing 


















FIG. 6. (Color online) (a) Simulated 521 spectra of e-mode 
for empty cavity (black solid lines) and YIG loaded cavity at 
Hsxt=0 (blue dashed lines). The hybridization of the cavity 
and FMR modes appears at /roJIea;t=412 mT (red solid lines). 
The coupling strength g as a function of (b) at VF = 85 
mm and £s = 15, (c) at 4 = 0.1 mm and Ss = 15, and 
(c) Sj at IT = 85 mm and 4 = 0.1 mm, with full numerical 
calculation (solid circles) and analytical results (solid lines) 
using Eq. (1121) (red solid lines). 
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FIG. 7. (Color online) 52i spectra measured as a function of 
external magnetic field and microwave frequency for a YIG 
(d=0.5 mm) loaded ID cavity. Dashed lines show the two 
strongest coupled magnon modes: the right one is coupled 
with /i-modes and the left one with e-modes. 


III. EXPERIMENTAL RESULTS AND 
DISCUSSION 


higher order w — ujcr terms with lo — ujcr = ^iocR based 
on the first order approximation, the transmission coef¬ 
ficient of the CMP for e-modes is found to be 


521 = 


1 + 


^^'CR 


i{u) i(^ljJ — U}pMR) — A.UJFMR - 

( 12 ) 


with 

‘^'cR — ^CR — — 1 + (Wc/wcr)^], 

Auj'cr = Aujcr = 

5^^ = (13) 

Similar to that for the /i—mode, the coupling strength 
g' in Eq. (1131) is verified by comparing its dependence on 
Is, I and Es with the results from the full-wave simula¬ 
tion as shown in Fig. HDb-d), respectively, where W =85 
mm and e,j=15 in Fig. |6Kb), Zs=0.1 mm and es=15 for 
Fig. He), and Is = 0.1 mm and TF = 85 mm for Fig. 
Hd). Without any adjustable parameters, the agreement 
between the full-wave simulation and the analytical ap¬ 
proximation is good. Although g' oc Z“F2 fgj- e—mode 
is the same as g oc fQj. mode, g' is linearly 

proportional to IJ and Es, differing from g. 


The CMP experiments have been performed by insert¬ 
ing an YIG sample into the proposed ID microwave cav¬ 
ity. Due to the difficulty in obtaining YIG samples with 
an exact diameter of 1.61 cm to satisfy the requirement 
of our ID model, a set of rectangular YIG samples with 
a dimension of about 5.5x 3.5 mm^ was used instead. In 
the first experiment, the YIG sample (with a thickness 
of d=0.5 mm) was inserted into the cavity with a total 
length of IT = 2l+ls = 85 mm and ^21 spectra were mea¬ 
sured by successively changing the magnetic bias iText 
from 300 mT to 600 mT with a step size of 0.1 mT. Fig¬ 
ure [3 maps the amplitude of ^21 as a function of the 
frequency and the magnetic bias. It has been found that 
a number of magnon modes in the YIG sample couple 
to the cavity. The multiple resonance structure could be 
attributed to the magnetostatic modes^i^ in our YIG 
sample, which have also been observed in previous cou¬ 
pled magnon-photon systemsi^d^. A carefully inspection 
indicates that there are two magnon modes (guided by 
the dashed lines) strongly coupled with the cavity modes: 
one coupled with /i-mode and another coupled with e- 
mode. The separation between them is only 6.0 mT. 

To simplify the complicated 3D wave propagation 
problem caused by (1) the lateral size (~0.2 cm^) of the 
YIG sample, which is much smaller than the cross section 
(~2 cm^) of the circular waveguide, and (2) spatially non- 
uniform precession of the magnon system, a filling factor 
^ = Is/d that converts the thickness d of the YIG sam¬ 
ple to be an effective Is in our ID model was introduced. 
As verified in the following experiments, it doesn’t af- 
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FIG. 8. (Color online) (a) S 21 spectra for /i-mode at 12.157 
GHz measured as a function of external magnetic field and 
microwave frequency. The dotted line refers to the pre¬ 
dicted FMR dispersion according to Kittel’s formula, while 
the dashed line is the cavity mode, (b) Typical 521 spectra for 
the CMP. (c) The calculated evolution of hybridized FMR and 
cavity mode corresponded to g/^n = 78 MHz with ^=0.18. 
(d) Dispersion of the hybridized CMP for measurements (sym¬ 
bols) and calculation (solid line) based on Eq. (HOI) . 
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FIG. 9. (Color online) (a) For h—mode S 21 amplitude map¬ 
ping as a function of the microwave frequency uj and the mag¬ 
netic bias for three different cases: (a)d = 0.5 mm and 
IF = 85 mm, (b) d = 1 mm and IF = 85 mm, and (c) d — 0.5 
mm and IF = 242 mm. (d) Typical spectra of S 21 shows the 
splitting of the photon-like and magnon-like CMP when they 
have equal amplitude. The deduced coupling strength are 78 
MHz, 121 MHz, and 48 MHz. 


feet the dependence of the coupling strength on either I 
or Is and furthermore it can quantitatively explain the 
experimental observation. 

The coupling of the magnon-photon system near the 
/i-mode at ujcrI‘2'^ = 12.157 GHz as shown in Fig. UKa). 
When the magnon modes are far from the cavity modes, 
they follow the dotted line, which is the dispersion cal¬ 
culated from Kittehs formula. However, when the pred¬ 
icated magnon frequency matches with the cavity res¬ 
onance frequency indicated by the dashed line, the cou¬ 
pling between microwave and magnetization dynamic hy¬ 
bridizes the magnon and cavity modes and a pronounced 
splitting of cavity modes is observed. 

Typical 521 spectra in Fig. [8Kb) illustrates the evo¬ 
lution of the hybrid magnon and cavity mode. When 
the resonance frequency difference between the magnon 
mode and the cavity mode is large, the magnon is only 
weakly excited by the electromagnetic field, while the 
amplitude of the cavity mode is orders of magnitude 
larger. The amplitude of magnon modes is significantly 
enhanced when it approaches the cavity mode. This 
effect is due to the fact that the coupling of magnon 
and cavity modes generates hybridized cavity photon- 
magnon quasi-particles, i.e., cavity magnon polaritona^^. 
At noHext = 381 mT corresponding to ujfmr = ^CR, 
the photon-like and magnon-like CMP have equal am¬ 
plitude and occur at ojcr ± 5, indicating that the mi¬ 
crowave energy is effectively transferred from the pho¬ 
ton system to the magnon system through the mutual 


electrodynamic coupling between them. Similar to the 
optical systemi^^ and metamaterial systent^S, the cou¬ 
pling between magnon and photon could experience elec¬ 
tromagnetic induced transparency and other coupling 
phenomenal^. 

For comparison, we numerically calculate the |52i| co¬ 
efficient using Eq. (flUl) and plot it in Fig. jSKc), which re¬ 
produces the measured coupling feature. The resonance 
frequencies of the hybrid modes are plotted as a function 
of the magnetic bias H^xt in FigjSKd), where the calcu¬ 
lated dispersion (solid line) agrees very well with the mea¬ 
sured data (symbols). To match g/2'K =78 MHz deduced 
from the measured separation between the the hybrid 
modes at g.oHext =381 mT, a filling factor of ^=0.18 is 
introduced to simplify the complicated 3D problem while 
other parameters are determined from independent ex¬ 
periments. 

To verify the scaling effect of Ig and I on the coupling 
strength g for the ft,-mode, the second experiment has 
been performed near 12.157 GHz by varying Ig and 1. Fig. 
|9Kb) presents the amplitude mapping of ^21 as a function 
of the frequency and the magnetic bias in the vicinity of 
/i-modes at 12.157 GHz, where a d=l mm thick YIG 
is loaded while keeping the cavity length IF = 85 mm 
unchanged. Similar dispersion of the CMP is observed, 
however, with a much larger separation compared with 
the dispersion measured for the YIG sample with (i=0.5 
mm. The enhancement of the coupling strength can be 
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FIG. 11. (Color online) Coupling strength as a fnnction of 
y/ls/l and Is^/ls/l for h-mode and e-mode, respectively. Sym¬ 
bols are experimental data with Is = and solid lines are 
calculation based on the analytic solution using fj,oHext = 380 
mT, ojcrI^tt = 12.2 GHz for the /i-mode and fJ-oHext = 420 
mT, cjch/Stt = 13.2 GHz for the e-mode, respectively. 


FIG. 10. (Color online)(a) For e—mode S 21 amplitude map¬ 
ping as a function of the microwave frequency uj and the mag¬ 
netic bias Hext for three different cases: (a) d = 0.5 mm and 
IF = 85 mm, (b) d = 1 mm and IF = 85 mm, and (c) d = 0.5 
mm and IF = 242 mm. (d) Typical spectra of 5'2i shows the 
splitting of the photon-like and magnon-like CMP when they 
have equal amplitude. The deduced coupling strengths g' are 
49 MHz, 130 MHz, and 24 MHz. 


clearly revealed by inspecting the 5'2i spectra at ujfmr = 
LocR as shown in Fig. [Sjb), where g/27r=121 MHz is 
deduced. When compared with the value of g/27r= 78 
MHz measured for a 0.5 mm thick sample, the obtained 
scaling factor is 1.5 which is close to the expected scaling 
factor of v^. 

We also checked the I dependence by elongating the 
length of the circular waveguide to W =242 mm but 
keeping c?=0.5 mm unchanged. As expected a weaker 
coupling was found as shown in Fig. [n(c). The mea¬ 
sured g/27r=44 MHz is close to the expected value of 
(7/271=48 MHz scaled by Therefore, the scaling 

effect of the sample thickness and the length of the cav¬ 
ity is demonstrated for h-modes by this experiment. We 
note that in this experiment the resonance frequency of 
the magnon mode slightly shifts with the sample thick¬ 
ness and the cavity also moves with 1. However, a care¬ 
ful inspection indicates that the resultant changes in the 
coupling strength should be less than 5%. 

In the third experiment, the scaling effect of Is and I 
on the coupling strength g' for the e-mode was verified. 
For the measurement performed for d = 0.5 mm and 
IF = 85 mm, the coupling strength is about 49 MHz for 
the e—mode at 13.20 GHz, which is weaker than that of 
78 MHz for the h—mode at 12.157 GHz. Increasing d to 
1 mm the coupling strength g' is raised to 130 MHz as 
shown in Fig. fTUl b). The estimated scaling factor of 
for e—mode is in agreement with Is predicted by our 


theory Eq. m- Similarly to the h-modes, the coupling 
strength g' scales approximately with using the val¬ 
ues of 49 MHz and 24 MHz for IF =85 mm and IF =242 
mm, respectively. 

To summarize our experimental observation and the¬ 
oretical expectation. Figure [11] highlights different de¬ 
pendence of the coup ling strength, where g is linearly 
proportional to \/ls/l while g' is linearly proportional to 
hy/ls/l- To fit the analytic results (solid lines) calcu¬ 
lated based on Eq. (HU and Eq. dn for the h- and 
e-modes, respectively, the filling factors for experimental 
data (symbols) are found to be ^ =0.18 for the h-mode 
occurring at ^ 12.2 GHz and ^ =0.32 for the e-mode oc¬ 
curring at ^ 13.2 GHz. The difference in ^ is believed to 
be due to the detailed wave propagation of the h-mode 
and e-mode. 


IV. CONCLUSIONS 

In conclusion, we have developed a transfer matrix 
approach to investigate the dispersion of the cavity- 
magnon-polariton (CMP), created by inserting a low 
damping magnetic insulator into a high-quality ID mi¬ 
crowave cavity. In the limit of long wavelength an an¬ 
alytical approximation of the CMP dispersion has been 
obtained, where the coupling strength characterizing the 
interplay between the magnetization and microwave dy¬ 
namics is determined by the sample thickness, the per¬ 
mittivity, the cavity length as well the the parity of the 
cavity modes. These scaling effects of the cavity and ma¬ 
terial parameters is then confirmed by experimental data. 
The analytic solution indicates the universality of the ID 
CMP, in which the coupling strength is not sensitive to 
the detailed design of the cavity. This work not only pro¬ 
vides a detailed physical understanding of the coupled 
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dynamic system that could be used to quantify the char¬ 
acteristic phenomena associated with different coupling 
regimes, but also could be potentially used for designing 
dynamic filters and switch devices for microwave appli¬ 
cations. 
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Appendix A: Details of the wave propagation in the 
cavity 


In this Appendix, we present some details of the ana¬ 
lytical solution for the wave propagation in the microwave 
resonance cavity, which can be used for the quantitative 
analysis of the coupling between FMR and cavity reso¬ 
nance. We start with the transfer matrix (the ABCD ma¬ 
trix in microwave terminology) following the con¬ 

vention. For a cascaded system of many components such 
as in FiglTJb), the transfer matrix between ports A and B 
is determined by multiplying the ABCD matrix for each 
component as 2 ^ 

(A1) 

The transfer matrix (Mi) of the air-filled waveguide 
section between port A and port L is 

M - ( cos(fcoO -iZosin(fco/)\ , . 

^ sin(fcoO cos{kol) J ’ ^ ' 

the transfer matrix (Mg) for the sample-loaded waveg¬ 
uide (yellow area between port L and R in Fig[TJb)) is 



section between port R and port B is 


M — I cos(kol) -iZosm(kol)\ 

^ (—iZg^sin(kol) cos(kol) )' 


(A4) 


In the above equations the wave number and the char¬ 
acteristic impedance are defined as 


ks 



and Zg = 


h 


(A5) 


Here ui, ujc, and £ 0 ( 1 * 0 ) are the frequency of the 

operating microwave, the cut-off frequency of the circular 
waveguide, the relative complex permittivity (permeabil¬ 
ity), and vacuum permittivity/permeability, respectively, 
c = l/y/Eofio is the vacuum speed of light. For empty 
waveguide (air-filled waveguide), the wave number and 
characteristic impedance are fco and Zq using Eg = fjg = 1. 

The scattering equations between the ports A and B 
are expressed as 


fS aa Sab\ ( e(4\ 
Sbb) V®b/ 


(A 6 ) 


Mg = 


cos{kglg) —iZg sin(A^Zs) 
—iZ^^ sin(A;j4) cos{kglg) 


(A3) 


and the transfer matrix (M 2 ) of the air-filled waveguide 

I 


where the where the superscript “—’’and “-I-”indicate 
that the direction of wave propagation is from a circular 
waveguide to the transition and vice versa. Thus the S- 
parameter between port A and B can be calculated as2£ 


Saa Sab\ 
Sba SbbJ 


/ A+BjZo-CZo-D 
( A+B/Zo+CZo+D 

\A+B/Zo+CZo+D 


2{AD-BC) ' 
/ Zq-\-C Zq-\-D 
-A+BIZq-CZq+D 
A-\-B/Zq-\-CZq-\-D , 


^\n{ksls){Zl-Zl) 


‘HZq Zs 


_ 2ikol I (ZJ + ZJ) sin{ksls)+^'iZoZs cos{ksls) (Zj+Zj) sin{ksls)-\-2iZQZs cos{ksls) 


2iZo Zs 


sin{ksls){Zi-Z^) 


(Zq+Z^) sin{ksls)-\-2iZoZs cos{ksls) (Zq+Z^) sin{ksls)-\-‘2iZoZs cos{ksls. 


(A7) 


Notice that the co-polarized TEll mode exists in 
the circular waveguide and the mono-polarized TEIO 
mode exists in the rectangular waveguide. The resulting 
boundary conditions for port A and B (in x — y coordi¬ 
nates) are 





(A 8 a) 

(A 8 b) 


where and are the incident waves from 

port A and B with x and y polarization, respectively. 
In the apparatus shown in Fig[T] y— and y' — polarized 
TEIl waves entering the transitions from port A and B, 
are reflected with a reflectivity R and phase jump of ipy, 


while X— and x'— polarized TEll waves are transmit¬ 
ted through the transition without reflection. Thus the 
reflecting matrix F aa and Tbb can be expressed as 

= (A9a) 

_ /cos( 0 ) — sin( 0 )\ / cos{0) sin( 0 )\ 

~ \^sin( 0 ) cos( 0 ) J \^— sin( 0 ) cos( 6 i )J ' 

(A9b) 

Combining the scattering equations Eqs. (IA 6 II between 
port A and B and using the boundary conditions Eqs. 
(IA 8 bl) at port A and B, it is found^^ 

A — Saa^aa —Sab^bb \ (_ fS aa Sab\ ( 
\ —Sba^aa ^ — Sbb^bb) Ve))/ \Sba Sbb) 

(AlO) 
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Taking into account the x — y co-polarized TEll wave, 
Eq. (lAlOt) describes the relation between the transmit- 

I 


ting/reflecting wave and incident wave in the x — y coor¬ 
dinate, which can be written as 


^Ay 

^Bx 

\^By) 



0\ _ (Saa 


0 

'y 0 Saa 
Sba 0 
0 Sba 




^AA 


_ fS ab 0 
0 Sab 
Sbb 




0 

Sbb 


bb 


(Saa 

0 

Sab 

0 

/Anc\ 
I^Ax \ 

0 

Saa 

0 

5ab 

Anc 

^Ay 

Sba 

0 

Sbb 


Anc 

^Bx 

\ 0 

Sba 

0 

Sbb) 

VByJ 


(All) 


As the experimentally accessible parameters are the By solving Eq. (lAllIl the final transmission coefficients 
transmission coefficients S '12 and/or 521 between ports 1 can be obtained: 
and 2 of the assembled waveguide cavity, defined as 


521 — 


-'Bx' 

„inc 


Cbx cos{0)-egy sin( 6 >) 


e 


me 

A 


(A12a) 


5 


^Ax 
^inc ■ 
-'B 


(A12b) 

I 


cos(g)5AB[-l + R^e^^^ySABSBA + e^^yRSeB - e^^yRSAAj-l + e^^^RSBB)] 
— 1 -I- cos^[9 )BSSabSbA + e’^’^yRSuB — e*‘^«i?5AA(—1 + e'-^^yRSuB) 


(A13) 


Appendix B: Wave propagation in the vicinity of the 
ferromagnetic resonance 

In this Appendix, we present some details of the ana¬ 
lytical solution for the wave propagation near the FMR. 
In the vicinity of the FMR frequency, the wave propa¬ 
gation is related to the dipolar field associated with the 
interaction of spins from Maxwell’s equations, written as 


V X h = EsEo 

V X e = —/xq 


ere. 


de 

m 

a(h-t-m 


dt 


(Bla) 

(Bib) 


where a is the conductivity of the ferromagnetic sample. 
For the insulating ferrite YIG a tends to be zero. 

Eliminating e from Eq. (IBla|) by taking the curl of it 
and substituting by Eq. (IBlbl) . we obtain 

V^h — EsEoyoAh = V(V • h) -I- EsEoyoAm. (B2) 


The motion of M is governed by the phenomenological 
Landau-Lifshitz-Gilbert equation^. 


dM 

dt 


- 7 (M X ha -k — 



(B3) 


Here 7 = yoe/rrie is the effective gyromagnetic ratio 
of an electron (charge e and mass me), a is the Gilbert 
damping parameter, Ms is the saturation magnetization. 
Inside the ferrite sample, the magnetization M = Mo-km 


and magnetic field = He^t+h consist of static compo¬ 
nents Mo and Hex* as well as the dynamic components 
m and h including a time dependent term of As¬ 

suming idext along y direction, from the vector relation 
of Eq. (IB3|) we can obtain the well-known Polder tensor 
X expressed as 


where 


f XL 0 
X= 0 0 

\-iXT 0 



(B4) 


XL 

XT 


lilHext - iauj)Mo 
(jHext - iaujY - ’ 

xMquj 

ixHext - iaujy - 


(B5a) 

(B5b) 


Using the approximation of long wavelength, i.e., the 
wave propagates only along the z direction, the cou¬ 
pling equation of Maxwell’s equations (IB2I) and Landau- 
Lifshitz-Gilbert equation (IB3I) become 

/esEqUouj'^^I + Xl) - iEsEoyoUj'^XT 

-lEsEoyoUj'^XT EsEoyoUj‘^{l + XL)J 

(B 6 ) 

The resultant wave number is 


kg — ; (BI’) 

where the effective permeability yeff is defined by yteff = 
[{l + XLr-xl]/{l + XL). 
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